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O . Abstract. Submillimetre (sub-mm) galaxies have very high infrared (IR) 

' luminosities and are thousands of times more numerous at z ~ 2 than local 

Q \ ultra-luminous IR galaxies. They therefore represent a key phase in galaxy evo- 

lution which can be missed in optical surveys. Determining their contribution 
to the global star formation rate requires dissecting their IR emission into con- 
\ tributions from starbursts (SB) and active galactic nuclei (AGN). There are 

■ several examples of AGN systems which masquerade as SBs in either the IR or 

' X-ray, and SBs can often look like AGN in some wavebands. A combination 

of SB and AGN emission is not unreasonable, given models of merger-driven 
evolution. To assess in detail what powers the intense IR luminosity of sub-mm 
04 ' galaxies it is important to obtain a complete multi-wavelength picture. Mid-IR 

spectroscopy is a particularly good probe of where the intense IR luminosity is 
coming from. We present the first results from a program to obtain Spitzer IRS 
spectroscopy of a sample of high redshift galaxies in the GOODS-N field, a large 
fraction of which are sub-mm galaxies. This field is already home to the deepest 
O ^ X-ray, optical, IR and radio data. We piece together the sub-mm data with 

the Spitzer photometry and IRS spectra to provide a well sampled IR spectral 
energy distribution (SED) of sub-mm galaxies and determine the contribution 
to the bolometric luminosity from the AGN and SB components. 
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1. Introduction 



Enormous effort has gone into understanding the nature of sub-mm selected 
galaxies since they were first discovered with the Submillimetre Common User 
Bolometer Array (SCUBA, Holland et al. 1999) just under a decade ago (e.g. Smail 
et al. 1997; Barger et al. 1998; Hughes et al. 1998). A general understanding of 
the sub-mm population and its role in galaxy evolution is severely limited by the 
poor spatial resolution of current sub-mm telescopes and the faintness of these 
sources at other wavelengths (e.g. Pope et al. 2005). The sub-mm population 
appears to consist of massive objects (Borys et al. 2005; Greve et al. 2005) at 
z ~ 2 — 3 (Chapman et al. 2005; Pope et al. 2005). These galaxies are rare and 
have very high star formation rates (SFRs, Lilly et al. 1999; Scott et al. 2002), 
and hence may represent an early phase in the evolution of massive elliptical 
galaxies. An outstanding question remains: are the incredible infrared lumi- 
nosities of sub-mm galaxies powered by huge bursts of star formation or AGN 
activity? To assess in detail what powers the intense IR luminosity of sub-mm 
galaxies it is important to obtain a complete multi-wavelength picture. 
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The sample of sub-mm galaxies discussed in this paper is taken from the 
GOODS-N SCUBA super-map (Borys et al. 2003; Pope et al. 2005). This is a 
purely sub-mm selected sample and it is almost completely identified at other 
wavelengths thanks to the deep multi-wavelength data available in the GOODS- 
N field (see Pope et al. 2006 for more details). 

2. Mid-infrared photometry 

The addition of the MIPS photometry to the multi-wavelength dataset for sub- 
mm galaxies provides powerful constraints on the shape of the SED, the nature 
of the power source, and the source redshift. In the left panel of Fig. 1, we 
plot the S24nm/ Ss5o^im flux ratio as a function of redshift as compared to several 
models. In general, this flux ratio as a function of redshift is lower than that of 
other ultra luminous infrared galaxies (ULIRGs) of the same luminosity at low 
and high redshift. This suggests that sub-mm galaxies either have higher levels 
of extinction or cooler dust temperatures. This can be tested further by looking 
in detail at the mid-IR spectrum. 

Using deep multi- wavelength follow-up observations of sub-mm galaxies, it 
is possible to study their SEDs to derive fundamental properties such as dust 
properties and SFRs. In the right panel of Fig. 1, we plot a composite rest- frame 
SED for the sub-mm sources with spectroscopic redshift as compared to several 
models. We find that sub-mm galaxies have cooler average dust temperatures 
than those of local ULIRGs of the same luminosity (see Pope et al. 2006 for more 
details). This result is confirmed with the addition of deep 70 /im imaging of 
these sub-mm galaxies (Huynh et al. 2006). This may indicate that their far-IR 
emission is more extended than that of local ULIRGs, in which the majority of 
the IR emission comes from within the central kpc (Charmandaris et al. 2002). 
The SEDs of sub-mm galaxies are also different from those of their high red- 
shift neighbours selected at near-IR wavelengths (e.g. BzK galaxies, Daddi et 
al. 2005), whose mid-IR to radio SEDs are more like those of local ULIRGs. We 
may understand this difference in terms of the evolutionary scenario advocated 
by Vega et al. (2005), in which enshrouded star forming galaxies undergo four 
distinct phases characterized by different mid- and far-IR colours. Under their 
scheme, the sub-mm galaxies, with cooler temperatures, may be an earlier phase 
in the star formation than BzK galaxies. 

3. Mid-infrared spectroscopy 

Many sub-mm galaxies show evidence of an AGN, as determined either through 
X-ray observations (Alexander et al. 2005), or with optical spectra (Swinbank et 
al. 2004; Chapman et al. 2005), however it is unclear if the AGN is a significant 
contributor to the bolometric luminosity of the galaxy. Mid-IR spectroscopy can 
help determine what powers sub-mm galaxies since these wavelengths contain a 
number of features which are sensitive to the nature of the energy source (Clavel 
et al. 2000; Hudgins & Allamandola 2004). 

We obtained Spitzer IRS observations of a sample of the sub-mm galaxies 
and AGN in GOODS-N (G02 PI: R. Chary). Fig. 2 shows the rest-frame raw 
IRS spectra for 10 sub-mm sources from this sample, offset for clarity. All of 
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Figure 1. (left) Evolution of the S^A^ml Ss50tim ratio with redshift for sub- 
mm sources in GOODS-N as compared to several models (see Pope et al. 2006 
for more details). The open circles show the average values of the sub-mm 
sources in four redshift bins, with roughly equal numbers of sources in each. 
The filled triangle is the average colour found for BzK galaxies (Daddi et 
al. 2005). The dotted curve is a cool SED template from Chary & Elbaz 
(2001, CE01 hereafter) models, the solid curve is the same CE01 template with 
additional additional extinction from Draine (2003), the dashed curve is the 
observed SED of Arp220 and the dash-dot curve is Mrk231, an IR-luminous 
AGN. The shaded region represents the redshifts where the 9.7 /um silicate 
feature passes through the 24 jim. passband. The generally lower S^/S^o 
ratios for the sub-mm sources suggest that they either have higher levels of 
extinction or cooler dust temperatures, (right) Composite rest-frame SED for 
GOODS-N sub-mm sources with spectroscopic redshifts. The solid curve is 
a modified CE01 model for a cool starburst galaxy, where we have applied 
additional extinction using the Draine (2003) models. The dash-dot curve is 
a CE01 model with the same luminosity but a warmer dust temperature. 



these spectra show at least weak polycyclic aromatic hydrocarbon (PAH) emis- 
sion. We have obtained several new redshifts from these IRS spectra (Pope et 
al. in preparation). Fig. 3 shows the mean IRS spectrum for these sub-mm 
galaxies compared to those of several local galaxies. This composite spectrum 
was calculated by taking the mean of the individual spectra. Note that using the 
median instead of the mean does not make a significant difference to the shape 
of the composite spectrum. We have excluded the top source in Fig. 2 since its 
steeply rising mid-IR continuum is not seen in the other 9 spectra in our sample. 
This source is not typical of our sample and is likely to have a stronger AGN 
component. In the composite spectrum, we see the 6.2, 7.7, 8.6 and 11.3 //m 
PAH features much more clearly. While Arp220 is often considered the typical 
local analog to high redshift sub-mm galaxies, it appears to have more silicate 
absorption and less 6.2 /im PAH emission than the average sub-mm galaxy. In- 
terestingly, the mid-IR spectra of the sub-mm galaxies seems to resemble that 
of M82 with the addition of a shallow power-law component (~z^~ ' 3 ). The 
total contribution of the power-law component to the total luminosity in this 
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Figure 2. Raw IRS spectra of sub-mm galaxies in GOODS-N. Each spec- 
trum has been offset for clarity and the horizontal lines show where zero is for 
each. Vertical lines indicate the positions of the PAH emission features and 
the 9.7 fxm absorption feature. 



wavelength range is ~ 30%. M82 is a typical local starburst galaxy but it is at 
least 2 orders of magnitude less luminous in the IR than the average sub-mm 
galaxy. This composite suggests that the mid-IR emission in sub-mm galaxies 
is dominated by star formation activity and there is little contribution from an 
AGN. The average 7.7/um line to continuum ratio for this sample is ~ 2.5-3.5. 
Comparing this with Figure 4 from Genzel et al. (1998), we find that this ratio 
for sub-mm galaxies is consistent with the ratios for SBs and ULIRGs and not 
as low as those found for AGN. Our preliminary analysis suggests that sub-mm 
galaxies are powered primarily by SBs and not AGN activity. 
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4. Future Work 

We are currently working on fitting templates to quantify the level of AGN con- 
tributing to the IR luminosity and investigating various emission and absorption 
line strengths as diagnostics for the star formation activity. Combined with the 
other multi-wavelength data that are already tabulated on these sub-mm sources, 
we will create a complete picture of the energetics of sub-mm galaxies (Pope et 
al. in preparation). These data will help determine if sub-mm galaxies and AGN 
populations are connected via an evolutionary sequence. 
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Figure 3. Composite IRS spectrum of sub-mm galaxies (Pope et al. in 
preparation). The dashed curve is M82 (Forster Schreiber et al. 2003) plus 
a shallow power-law component (^i/~ - 3 ) while the dotted curve is Arp220 
(Charmandaris et al. 1999; Fischer et al. 1999). The average sub-mm galaxy 
shows PAH features, some silicate absorption and a fairly flat continuum in- 
dicating little or no contribution from an AGN to the luminosity at these 
wavelengths. 



